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ABSTRACT 

Pseudo-Goldstone bosons couple with photons through a P and T violating inter- 
action of the form Cx = g a j-y aFF. Strong magnetic fields in rotating compact 
stars induce a non-zero E ■ B outside the stellar surface which acts as a source 
for the pseudo-scalar field. Pulsar signals propagating through this pseudo-scalar 
'hair' suffer a differential time lag between the left and the right circularly po- 
larised modes because of the P and T violating pseudoscalar boson-photon in- 
teraction. Determination of this time lag upto a microsecond accuracy can lead 
to a measurment of (or rule out) pseudo-Goldstone boson-photon coupling upto 
g ail < 4.6 x 10" 9 GeV' 1 . 
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Goldstone bosons arise in particle physics models when a global symmetry 
is spontaneously broken. These include axions [1] which arise in many different 
versions of the Peccei Quinn symmetry breaking and may be the solution of the 
strong CP problem. Other examples include Majorons [2] which arise when the 
lepton number is broken spontaneously and are essential in order to reconcile pos- 
sible heavy neutrinos with astrophysical bounds (a neutrino heavier than 25eV 
and lighter than a few GeV has to decay into lighter neutrinos and Majorons in 
order not to overdose the universe). Goldstone bosons also arise in superstring 
theories [3] in the same multiplet as the graviton and the dilaton and it is pos- 
sible that some of the above mentioned Goldstone bosons are string axions. The 
coupling of Goldstone bosons with other fields is determined by the energy scale 
of the global symmetry breaking. Generically psedo-scalar boson-photon coupling 
can be written as C ail = g ail aF fJjV F tJ ' v with g ail ~ C/v where C is model de- 
pendent constant and v is the global symmetry breaking scale. Strong magnetic 
fields outside compact rotating stars [4] induce a non-zero E ■ B density outside 
the stellar surface which acts as a source for the pseudo-scalar field. Pulsar sig- 
nals propagating through this pseudo-scalar 'hair' suffer a differential time lag 
between the left and right circularly polarised modes because of the Parity and 
Time-reversal violating pseudo-scalar photon interaction. We show that for a typ- 
ical millisecond pulsar the time lag is 5t ~ 10 _6 sec for a coupling of the order 
g ail ~ 10~ 9 GeV. Measurement of such a time lag between the two circularly 
polarised modes of the order of 10 _6 sec have been reported in literature [5,6,7]. 
Due to the large systematic error involved these numbers should be regarded as an 
upper bound on St which translates to an upper bound on the pseudo-Goldstone 
boson-photon coupling, g ail < 4.6 x lO^GeV -1 . Our bound is strictly valid for 
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massless pseudo-scalar fields. We have estimated that this bound continues to 
hold for pseudo-scalars with inverse mass greater than the radius of the neutron 
star i.e. for m a < 10~ 10 ey. Hadronic axions have a non-zero mass (induced by 
QCD instantons) in the range of 10 -5 — 10~ 3 eV, and our bound does not hold 
for them. Previous astrophysical bounds have considered the coupling of axions 
to matter in the cooling rates of the Sun, the red giant stars and the Supernova 
1987A [8,9,10]. The best astrophysical bound for the symmetry breaking scale 
as obtained from the cooling rate of red giant stars is v < lQ~ s GeV~ l . Terres- 
trial experiments have attempted to produce axions by laser beams in a cavity with 
magnetic field. The upper bound on axion-photon coupling from such experiments 
is g ail < 2.5 x 10~ 6 GeV~ 1 for m a < 10~ 4 eV [11]. Experimental searches for 
axions in vicinity of the Sun depend on the conversion of axions to X-rays in the 
stellar magnetic field. Based on a survey of the X-ray flux from the Sun, Lazarus 
et al [12] give an upper bound of g ail < 3.6 x lO~ 9 GeV~ l for m a < 0.03eF. 
Our bound for g ail is of same order as that of Lazarus et al. In pulsar signals the 
time lag between the two circularly polarised modes for the current upper bound 
g ail ~ 10~ 9 Gey _1 is St ~ 10 _6 sec. If it is possible to analyse the pulse profiles 
at a resolution higher than 10 _6 sec then it may be possible to actually observe a 
pulse splitting between the left and the right polarised modes and thereby deter- 
mine the value of g ail - The time lag also depends upon the frequency (St oc uj~ 2 ) 
and by analyzing pulsar signals at frequencies at lower than 100 MHz (for which 
St(100 MHz) = 10~ 6 sec) it may be possible to lower the bounds on g ail or make 
a positive measurement of pseudoscalar-photon coupling. 

Interaction of pseudo-Goldstone bosons with photons is described by the La- 
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grangian 

£ = \{d,a){d»a) - - \g ail aF, v F^ (1) 

where F^ u = (1/ \/ ri g)e^ upa 'F pa is the dual of F^ u . The equation of motion for 
axion field given by (1) is 

(-j=)d ll W=gd>>a) = -iSTLF^. (2) 

In terms of the electromagnetic field strengths E^ = F Q i and B % = e l i k Fjk, the 
source term in the right hand side of the equation (2) is g ail {l/ \/ ZI g){E ■ B). In 
order to have an non-trivial axion field configuration outside the compact stars 
we must have an external non-zero pseudoscalar density E ■ B. In dyonic and 
Kerr-Newmann family of blackholes [13] the external E ■ B is generated by the 
electric/magnetic charges and rotation and acts as a source for the pseudo-scalar 
hair in black holes [13]. Pulsars have a strong dipolar magnetic field with field 
strengths Bq ~ 10 12 Gauss which is determined by measuring the frequency of 
synchrotron X-ray emissions from their neighbourhood. In the oblique rotor model 
[14] for pulsar emission the dipole moment jl of a pulsar is not aligned along the 
rotation axis O (while in the case of blackhole they are parallel) . The dipole axis 
precesses around the spin axis with an angular frequency Q and radiates energy 
at the rate E = (-2/3) | 'jl |= -(1/6)S 2 R 6 Q A sin 2 a ~ 10 38 er# sec' 1 for a 
typical pulsar. The angle a can be determined by measuring O and O for a given 
pulsar. In the aligned rotor model [14] the angle a ~ (a small misalignment is 
necessary to produce the pulsation) and rate of radiation E ~ — S 2 i? 6 4 . The 
synchrotron radiation from the plasma accreting along the magnetic poles appears 
as the pulsed signal to an observer in the cone swept out by the radio beam as 
it precesses about the spin axis. The instantaneous magnetic field at point r is 
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given by 



B = ^m^)-v (3) 



where | p, |= BqR 3 and So is the maximum magnetic field strength at the surface 
(r = R) of the pulsar. Typically B ~ 10 8 - 10 12 Gauss and R ~ 12 km. Magnetic 
field in a rotating frame appears as an electric field 

E = -Ox fx B. (4) 

The magnetic field (and consequently the electric field) has a time dependent part 
oscillating with frequency O generated by the component of \x perpendicular 
to Q which precesses with angular frequency O. The time independent part of B 
and E are generated by the parallel component of \x. The time average of field 
strengths at time scale > O -1 is determined by < fl >= fj,\\ = B R 3 cosa (where a 
is the angle between jl and ft). The average value of electric field outside the stellar 
surface is then determined by (3) and (4) by matching the boundary conditions at 
r = R. Since E = in the highly conducting interior region, the component of E 
parallel to the surface is zero at the boundary and the perpendicular component of 

— * 

E changes discontinuously at r = R. The average electric field outside the stellar 
surface (r > R) is given by [15] 



B R 5 ttcosa 



2. 



< E ext (r) >= — , 3(sin 2 9 - -)f - 2sin9cos9 9 . (5) 



r 1 r v :>> 
The average value of the pseudoscalar density E ext ■ B is 

p8 

< E ext ■ B >= -B 2 ncosa -^cos 3 9 (6) 

r' 

which appears as the source term for the pseudo-scalar fields in equation (2). In 
the absence of pseudo-scalar coupling the strong electric field (5) which is of the 
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order E ~ 10 8 volts/cm (for O = lsec~ l and So = 10 12 Gauss) will create a 
space charge of particles with charge density p = (1/4/t) V • E = —(l/n) B ■ O 
= (B R 3 Q/2nr 3 )(sin 2 6 - 2cos 2 9). This space charge is of the order of 10 11 cm -3 
and could effectively screen the electric field. Because of the g ail aFF coupling 
however the electric field is screened by the pseudo-scalar field and not primarily 
by the ionised plasma. This follows from the equation of motion of electromagnetic 
fields modified by the pseudo-scalar coupling term 

DpF"" = 0a77 («9 M a)F^, (7) 

— * — * 

which in terms of E and B gives the modified Maxwells equations 

V.E = -g ail (Va)-B, (8a) 

-d E + V x B = g ail (Va) x E + (d a)B . (86) 

The remaining two Maxwells equations follow from the Bianchi identity d^F^ = 
which gives 

V- S = 0, (9a) 
d B + Vx£ = 0. (96) 
From equation (8a) it follows that the effective field experienced by the plasma is 

E eff = E + g ail aB. (10) 

— * 

Therefore a nonzero a and B will anti-screen the electric field. The space charge 
induced in the plasma will neutralise E e ff in (10) but even with the plasma po- 
larisation E can remain non-zero. The standard picture of pulsar magnetosphere 
[15] is no longer valid in the presence of pseudo-scalar fields. This may have other 
observable consequences besides the one we study here. 
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We assume the geometry outside the surface of the pulsar is described by the 
Kerr metric 

ds* = (i - ™M.)<n* - (i - ^)~V 2 - r\ d e* + ^W) 

sin 2 6d(j)dt. 11 

r 

The form (11) is an approximation to the Kerr metric, valid as long as the an- 
gular momentum per unit mass A « GM. For a typical millisecond pulsar 
(M = M Q7 R = 12 km), we have A = 10" 6 km and GM = 1.5km. Therefore 
(11) is an adequate approximation to the Kerr metric. For the metric given by 
(11) the equations of motion (2) for the axisymmetric, time independent solution 
a(r, 9) reduces to 



l_d_ 



r 1 1 



2GM\ d 



1 d 



d 



)Yr a ^ d) \ + ^^em[ smd m a ^ d) \ 

= g ail (B 2 R 8 Q coso) — cos s 9. (12) 

The Greens function which is solution of V 2 G(x,y) = — 5 3 (x — y)/^/—g for the 
Laplacian operator in (12) is given by [13] 

G(r, S, r , fl ) = ESo QP, - l)0, - 1)P, (A), r > r„ 

= ES » c ' ,Q '(5F- 1 ) p, (^- 1 ) p ' (A) ' r<r <" (13) 

where A = cos6cos6 + sin6sin6 , C\ = (21 — 1)/4tcM and P; and Qi are the two 
associated Legendre functions. The solution of the inhomogenous equation (12) is 

a(r,6) = - dr d6 #orS» 2 floG(r,ff,r ,^) 

JR JO JO 



x i^g ail B 2 R 8 cosa O j 7 u . (14) 



cos 3 # 



Substituting (13) in (14) and evaluating the integral we obtain the Goldstone 
boson configuration outside pulsar (r > R) 

^-(^) * m ? y~ g+°(?)- <*> 

Propagation of electromagnetic radiation through this coherent pseudo-Goldstone 
boson field is governed by the wave equation obtained by combining the Maxwells 
equations (8) and (9) 

V^VB = - <7 arY (Va) x d B. (16) 

With the eikonal ansatz 

B(x,t) = B e lS ^ x ^ (17a) 

and 

-id^B = (d^S)B = k^B (17b) 

the wave equation reduces to 

k^B = i g B77 (Va x ujB) (18) 

where /c M = (a;, k) is the photon four momentum and B is the magnetic field 
of the radiation. The pulsed radiation from pulsars originates from the polar 
regions of the pulsar [16]. In this region a ~ ao(l/r 2 ) and Va = — (2ao/r 3 ) f is 
radial. Consider right and left circularly polarised combinations for the transverse 
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magnetic field B 

£± = (^±^) (19) 

where Bi and Bj are the components of B along any two fixed directions and ij 
orthogonal to f. Substituting (19) in (18) we find that the equations for the two 
circularly polarised modes B± decouple as 

k^B± =f g ail {d r a) u B± = 0. (20) 
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This shows that the propagating modes in a background with non-zero pseu- 
doscalar gradients are circularly polarised modes [17] with different phase veloci- 

— * 

ties. The dispersion relation for the B± modes propagating radially from the poles 
is 

2 / 2GM\-i 2GM- 



(l-^— ) -fe^l-— — J = ±g ajJ (d r a)k . (21) 
The Eikonal phase can be written as 

S(r, t)=ut- J k r dr. (22) 

The propagation time t to reach a point r is given by the orbit equation dS / du = 
which gives 

d r 

t = — / k r dr. (23) 
Substituting for k r from (20) we have 

where £+(£_) are the arrival times of the +(— ) circularly polarised modes. Substi- 
tuting for d r a(r) from (15) (with cosO = 1) and evaluating the integrals in (24) we 
obtain the arrival times of the pulses from a pulsar at distance r from the earth 

(± = r + 2G M in (-) ± -(-) g A _ . (25 ) 

The second term is the familiar time delay due to the gravitational field of the 
pulsar and the last term shows the birefrigence properties of the pseudo-scalar 
background. The differential time delay between the positive and the negative 
polarised modes is 



5V575/ u 2 (GMf 
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Pulsar signals are recorded in the radio frequency range as low as ui ~ 100MHz. 
For a typical millisecond pulsar M = M Q = 10 57 GeV, B = 10 12 Gauss ~ 
10~ 9 GeV 2 , R = 12 km = 6 x 10 19 GeV _1 , cosa ~ 1 the differential time de- 
lay turns out to be 10 -6 sec for g ail ~ Vd~^GeV~ x . By putting an experimental 
bound on St one can therefore obtain the following upper bound on the pseudo- 
Goldstone boson-photon coupling. 

9 - = 46 x 10 GeV (tfw) h^)hr) (-S-) 

/ \V2, M ^3/2 

X llOOM^J Im©J • (27) 
Loseco et al [5] have analysed the signals from the millisecond pulsar PSR 
1937+214 using the data of Cordes and Stinebring [6]. They report that a time lag 
5 x 10 -6 sec is observed between the positive and the negative polarised modes. 
Klien and Thorsett [7] report a differential time lag of 0.37±0.67//sec. At the one 
sigma limit one can put an upper bound on the time lag to be St < 10 -6 sec. This 
leads to the upper bound on the goldstone boson photon coupling 

g ail < 4.6 x lO^GeV' 1 . (28) 

From (26) it is clear that by going for lower frequency pulses one should 
observe a larger time lag since St oc u~ 2 . Therefore by analysing lower frequency 
signals it should be possible to observe a non-zero time lag which would be a 
positive signature of pseudo-scalar photon interaction. 
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